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We estimate the production rate of axion-type particles in the core of the Earth, at a temperature 
T « 5000 K. We constrain thermal geo-axion emission by demanding a core-cooling rate less than 
0(100) K/Gyr, as suggested by geophysics. This yields a "quasi- vacuum" (unaffected by extreme 
stellar conditions) bound on the axion-electron fine structure constant a'^Y ^ 10~^*, stronger than 
the existing accelerator (vacuum) bound by 4 orders of magnitude. We consider the prospects for 
measuring the geo-axion flux through conversion into photons in a geoscope; such measurements can 
further constrain aSY . 



A variety of scenarios for physics beyond the Stan- 
dard Model (SM) give rise to light pseudo-scalar parti- 
cles, generically referred to as axions. The Peccei-Quinn 
(PQ) solution to the SM strong CP problem provided 
the initial context for axions . Axion-type particles are 
ubiquitous in string theory constructs and have also been 
considered in cosmological model building |^ . There are 
stringent astrophysical and cosmological constraints on 
the couplings of axions, as a result of which they are 
largely assumed to be very weakly interacting. Some of 
the strongest bounds on axion-SM couplings come from 
astrophysics, where stellar evolution and cooling argu- 
ments imply that the axion (PQ) scale / ^ 10^ GeV. 
Such analyzes are based on the requirement that new 
exotic processes should not significantly perturb a stan- 
dard picture of the energetics that govern the evolution 
of various astrophysical objects. Since axions (or other 
light weakly interacting particles) can directly drain en- 
ergy out of such objects, one can obtain bounds on the 
coupling of axions to matter. For a concise summary of 
various astrophysical bounds, see Ref. More recent 
astrophysical bounds on axion-type particles have been 
presented in Ref. Q. 

In this work, we consider the possibility that the hot 
core of the Earth can convert some of its thermal energy 
into a flux of axions of 0{eV) energy Then, it would 
be interesting to find out what bounds can be obtained 
from geological considerations and also to determine the 
prospects for discovering the geo-axions emanating from 
the terrestrial core. The Earth's core is at a temperature 
of around 5000 K corresponding to 0.4 eV. Although this 
is a much lower temperature than those of stellar interi- 
ors, which have temperatures of order keV, there are a 
number of considerations that motivate our analysis. 
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^ Non-thermal geo-axions, produced in radioactive decays within 
the Earth, have been examined in Ref. B]. This work does not 
find a currently detectable signal even in the most favorable case 
considered therein. 



First of all, the core of the Earth is only a short dis- 
tance away, compared to any astronomical object. This 
greatly enhances the prospects for measuring a geo-axion 
flux and can potentially compensate for the low core tem- 
perature. Secondly, the Earth's core is quite different 
from other axion emitting environments, being mainly 
made up of hot molten or crystallized iron. Hence, in 
principle, the intuition and calculations that apply to 
stellar plasmas may not be adequate to estimate geo- 
axion emission and new effects may need to be consid- 
ered. Finally, the Earth's center is a far less extreme 
medium compared to stellar media. The possibility of 
the dependence of axion properties on the environment 
has been proposed in the context of reports of large 
vacuum birefringence by PVLAS This result if it 
were confirmed would have implied an axion like parti- 
cle with an axion-photon coupling in stark violation of 
astrophysical bounds. As a consequence, a number of 
models were developed to reconcile the laboratory re- 
sult with the astrophysical bounds ^ . Although the 



initial PVLAS result could not be reproduced |I0 



highlighted the necessity for obtaining complementary 
bounds on axions in a wide variety of production environ- 
ments. If axion couplings are temperature and/or density 
dependent, the geo-axion bounds could be viewed as in- 
dependent new data on axion physics in "quasi-vacuum" 
conditions. Thus, in this letter we do not attempt to 
supersede existing, stringent astrophysical bounds, but 
to supplement them by examining axion production in a 
novel environment. 

Motivated by the above discussion, we will next derive 
an estimate for the thermal geo-axion flux. We will use 
geodynamical considerations to constrain this flux and 
hence the axion-electron coupling aae in the core. This 
bound is not competitive with its astrophysical counter- 
parts, but, as mentioned before, is derived in a very dif- 
ferent regime. Note that collider bounds on aae that are 
derived in a similar regime are much weaker than our 
geo-axion bound. We will then consider detection of the 
geo-axion flux, via magnetic conversion into photons, us- 
ing a "geoscope," in analogy with the helioscope concept 
Jill , . A discussion and a summary of our results are 



presented at the end of this work. 

The core of the Earth is mainly made of iron (Fc). The 
inner core, which extends to a radius of Ric w 1200 km, 
is thought to be in sohd crystaUine form at a tempera- 
ture T 6000 K. The outer core, which extends to Rc ~ 
3500 km, is made up of molten iron at T - 4000 K ^S). 
Since Fe is a transition metal, with the electronic con- 
figuration [Ar] 4s^, both 3d and 4s electrons are 
important in determining its properties. However, for 
a simplified treatment, we only consider the 4s electrons 
as nearly free. The effective nuclear charge seen by the 
4s electrons is Z^s ~ 5.4 @. 

Given that the solid iron core makes up a negligible 
mass of the total core, we will ignore its contributions to 
our estimate. This is partly done to avoid a complicated 
treatment of the interactions of electrons and phonons in- 
side a hot crystal, far from the plasma regime. However, 
we note that a more complete analysis should take these 
effects into account. We adopt « 5000 K« 0.4 eV |l| 
as the mean temperature of the molten iron core. We 
are also ignoring the contribution of other trace elements, 
such as nickel, which have more or less the same proper- 
ties as iron, for our purposes. Given the metallic nature 
of the core, we will treat it as a plasma composed of 
a degenerate gas of free electrons, with a Fermi energy 
Ep ~ 10.3 eV [ p^ . Th e resulting Fermi momentum is 
given hy pp = v'2™e Ep « 3.3 keV, where me ~ 0.5 MeV 
is the mass of the electron. These free electrons move in 
the background of Fe ions with effective charge Z^s — 5.4. 
The free electron density in the core is given by Pq] 



37r2 



(1) 



and hence we get rtg « 2 x 10^'^ cm ^. Let us define the 
radius 



/ 3 



V 47m£ 



1/3 



(2) 



for the mobile charged particles in the plasma, which we 
take to be electrons here. The quantity 



Z? 



apTn 



(3) 



with a = 1/137, is a measure of the relative strength 
of Coulomb interactions and the kinetic energy of the 



electrons. For the core parameters, we get a,. 



10- 



and r ~ 10^. We take F ^ 1 as indicative of a strongly 
coupled plasma [ p^ . Since the iron core of the Earth is in 
a molten state and not yet a crystal, this interpretation is 
reasonable, despite the large value of F. The effect of the 
geomagnetic field in the core on the density of states close 
to the Fermi surface can be neglected since the thermal 
energy is large compared to the energy difference between 
successive Landau levels. In any case, we note that a 
more detailed numerical treatment may reveal important 
corrections to the estimates that follow. 



Interestingly enough, there is an astrophysical environ- 
ment that is described by the above key features. This is 
the interior of White Dwarfs (WD's) which is a strongly 
coupled plasma of Carbon and Oxygen, supported by a 
degenerate gas of electrons, similar to the iron core of 
the Earth. Hence, we adopt the formalism used for WD 
cooling by axion emission in the bremsstrahlung process 
eN{Z,A) eN{Z,A)a |l|l, in order to estimate the 
geo-axion flux; Z is the ionic charge and A is the atomic 
mass. We will ignore Primakoff ||T^ contributions to this 
flux, resulting from the interactions of thermal photons 
in the plasma. This is justified, since the density of such 
photons is roughly given by (eV)'^ ~ 10^*^ cm^'^, which is 
much smaller than rig in the core. 

For a plasma with only one species of nuclei, the energy 
emission rate, in axions, per unit mass is given by p7[ 



{Z'a'aae)/{Amimu)T^i{pp), 



(4) 



where iriu — 1.7 x 10"^* g is the atomic mass unit and 
^{pf) is a numerical factor which only depends on pp. 
Numerical calculations relevant for WD's indicate that 
^ ~ 1 to a good approximation, over a wide range of 
parameters in the strongly coupled regime [p^ . We thus 
take ^ '--^ 1 in our calculations. For geo-axion emission, 
we then obtain 



Ea ^ lO"^ aaeT^ Crg g ^S" 



(5) 



where we have set Z — Z^g — 5.4, A — 56, and T3 = 
T/10'^ K. Given a core mass density of pc — 10 g cm~'^ 
and Tc ~ 5T3, we get 



$a ^ lO^^Qfae CrgS" 



(6) 



for the flux of geo-axions. 

It is interesting to inquire how geological considera- 
tions can constrain the estimate in Eq. (ph. As a simple 
criterion, and in the spirit of analogous considerations 
for stellar objects, we will demand that the rate of core- 
cooling $a be less than that inferred from geodynamical 
considerations. This rate has been estimated to be in 
the range of 100 K/Gyr= lO"'^ K/yr ||l§. Given that 
the heat capacity of the Earth's core is estimated to be 
~ 10'^'* erg/K ||l^, we get for the geological rate of 
core-cooling 



10 erg/yr '--^ 10 ergs 



(7) 



in agreement with Ref. ||l3[. Requiring < yields 

a® < 10~^* (core-cooling). (8) 

This bound is not strong compared to those from as- 
trophysics. For example, the bound from solar age is 
ctae ^ 10~22 and the one from red giant constraints is 
ctae ^ 10~2^ |l6|. However, the bound in (||) is within a 
few orders of magnitude of the solar. Again, we note that 
the bound in (H) is valid for a quasi-vacuum regime and 
not the extreme stellar environments. The closest such 
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bounds for quasi- vacuum environments are from e+e" 
collider experiments, and correspond to aae ^ 10^^^ 
(/ > 1000 GeV) weaker than our geo-axion bound 
yQ) by 4 orders of magnitude. Next, wc will examine the 
prospects for detecting a geo-axion flux consistent with 
this bound. 

The geo-axion flux Fa, corresponding to (f>a in Eq. (|^), 
at i?® ~ 6.4 X 10^ km (surface of the Earth) is given by 

Fa = $a/(4^ii'e) aaelO'° eVcm^^s-i. (9) 

Assuming average axion energy (Ea) — 1 eV, we get 

'^^'^ - aaelO'" cm-^s-i 



dAdt 



(10) 



for the flux of eV-axions at the surface of the Earth^ . 

In principle, there are two ways to detect these ax- 
ions: The obvious choice is to exploit their coupling to 
electrons aea via the axio-electric effect, in analogy with 
the photo-electric effect. The cross section for the axio- 
electric effect is reduced by a factor of [pOl 



Za \ nif. 



(11) 



compared to the ordinary photo-electric cross section. 
Using our bound derived in Eq. (^ and the typical size 
of photo-electric cross sections of 0(10"^") cm^, the re- 
sulting axio-electric cross section is 0(10"''^) cm^, which 
is quite small. For comparison typical neutrino cross sec- 
tions are O{10^'^^) cm . Therefore, we next consider the 
possibility to use magnetic axion-photon mixing, to con- 
vert axions into photons. We use the simple formula pM 



Pa 



BgajL 



(12) 



where B is the strength of the transverse magnetic field 
along the axion path, ga^y is the axion-photon coupling, 
and L is the length of the magnetic region. The above 
equation is valid when the q <^1/ L with 



q=\ml-wl\/{2Ea), 



(13) 



where ttIq is the axion mass and is the photon effective 
mass. 

Note that for ma ^ lO^'^ eV, the geo-axion oscillation 
length ^ 0.6 m in vacuum with = 0. We will 
assume this mass range for the purposes of our discus- 
sion. An obvious choice would be to consider an LHC- 
class magnet, such as the one used by the CAST experi- 
ment with _B = 10 T and L = 10 m. However, this 
magnet has a cross sectional area of order 14 cm^. Given 



^ For axions, the non-radial flux also contributes, thus in principle 
Gauf3' law can not be used. In this case, the difference compared 
to an exact treatment is about 5%. 



that the core of the Earth subtends an angle of order 
30° as viewed from its surface, we see that the CAST 
magnet will capture a very small portion of the relevant 
"field of view." Thus, we have to consider other magnets 
of similar strength, but larger field of view. Fortunately, 
such magnets are used in Magnetic Resonance Imaging 
(MRI), to carry out medical research. For example, the 
MRI machine at the University of Illinois at Chicago has 
a magnetic field of 9.4 T, over length scales of order 1 m 
[p2t . Hence, we will use B = 10 T and L = 1 m, as 
presently accessible values, for our estimates. 

The current best laboratory bound on (7^7 for nearly 
massless axions, derived under vacuum conditions, was 
recently obtained by the GammeV collaboration [p3| : 
Qa-y < 3.5 X lO-'^GeV"^ With this value the upper 
bound for the axion photon conversion probability is 



p 



3 X 10" 



(14) 



Thus, we conclude that magnetic detection is more 
promising than axio-electric detection. 
Using Eq. (|I^, we then get 



dNy 
dAdt 



(15) 



for the flux of converted photons in the signal. The bound 
in (^) then suggests that a sensitivity to a photon flux 
of order 1 cm^^s^^, using our reference geoscope param- 
eters, is required to go beyond the geodynamical con- 
straint and look for a signal. Modern superconducting 
transition edge bolometers have demonstrated single pho- 
ton counting in the near infrared with background rates 
as low as 10~^ Hz ||2^ and quantum efficiencies close to 
unity 1^ . With such a detector photon fluxes as small as 
10^^ cm^^ s^^ can be detected with an integration time 
of 10^ s. Our results on the prospects of direct search for 
geo-axions are summarized in figure |l|. 

Before closing, we would like to point out a few direc- 
tions for improving our estimates. First of all, our picture 
of the iron core is quite simplified. A more detailed treat- 
ment of electron-ion interactions in the molten core, as 
well as the inner core contribution, which was ignored 
here, could reveal extra enhancements or suppressions 
that were left out in our analysis. This could, in princi- 
ple, require a numerical simulation of the strongly cou- 
pled plasma (molten Fe) and the crystalline solid core. 
Another issue is the possible role that the Fe 3(i-orbital 
electrons play, given that they are delocalized over a few 
nuclei and may contribute to pseudo-scattering processes 
inside the hot Fe medium. Also, there could be impor- 
tant bound-bound and free-bound processes that result 
in the emission of axions from Fe atoms at high tem- 
peratures. These processes have been ignored here, but 
could provide contributions comparable to those we have 
estimated. In principle, more detailed geodynamical an- 
alyzes may yield stronger bounds on non-convective en- 
ergy transfer out of the Earth's core. This can result 
in tighter bounds on axion-electron coupling aae, in the 
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FIG. 1: The black lines show the resulting photon fluxes (p in 
units of cm~^ s^^ for a geoscope with L = 1 m and B = 10 T 
as a function of Uae and Qaj- The gray lines show the cur- 
rent astro-physical bounds on Oae from the cooling of white 
dwarfs Q and on Qa-y from the non-observation of solar axions 
by CAST The colored/shaded regions indicate the pa- 

rameter space excluded by photon regeneration and the 
cooling of the Earth (this work). 



regime we considered here. Finally, our estimate of a geo- 



scope signal assumed an axion flux transverse to the mag- 
netic field. Given the angular size of the core, as viewed 
through the geoscope, we expect that effective transverse 
field is, on average, suppressed by roughly cos^ 30°, which 
does not affect our conclusions, given the approximate 
nature of our estimates. 

In summary, we have derived estimates on possible 
emission of axions from the hot core of the Earth. Our 
analysis allows for possible dependence of axion proper- 
ties on non- vacuum production media, such as astrophys- 
ical environments. We approximated the molten core as 
a strongly coupled plasma of free degenerate electrons in 
the background of Fe nuclei. We adapted the existing ax- 
ion emission estimate from a White Dwarf interior, which 
is a strongly coupled plasma supported by a degenerate 
electron gas. We obtained the bound a®. < 10~^* on 
the axion-electron coupling, by considering geodynam- 
ical constraints on core-cooling rates. Given that geo- 
axions would originate from a far less extreme environ- 
ment than stellar cores, our bound is nearly a vacuum 
bound. Hence, our result improves existing accelerator 
constraints on aae, in vacuo, by 4 orders of magnitude. 
We also estimated the signal strength to be expected in a 
dedicated search for geo-axions using a geoscope, based 
on magnetic axion-photon conversion. 

Acknowledgments 

We would like to thank G. Khodaparast, S. King, and 
Y. Semertzidis for useful discussions. The work of H.D. 
is supported in part by the United States Department of 
Energy under Grant Contract DE-AC02-98CH10886. 



[1] R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38, 1440 

(1977); Phys. Rev. D 16, 1791 (1977). 
[2] J. Preskill, M. B. Wise and F. Wilczek, Phys. Lett. B 

120, 127 (1983); L. F. Abbott and P. Sikivie, Phys. Lett. 

B 120, 133 (1983); M. Dine and W. Fischler, Phys. Lett. 

B 120, 137 (1983); M. S. Turner, Phys. Rev. D 33, 889 

(1986) . 

[3] C. Amsler et al. [Particle Data Group], Phys. Lett. B 
667, 1 (2008). 

[4] P. Gondolo and G. Raffeh, arXiv:0807.2926 [astro-ph]. 
[5] Emission of axions from the core of Jupiter has been con- 
sidered in: C. Sivaram, Earth, Moon, and Planets 37, 155 

(1987) . 

[6] A. Liolios, Phys. Lett. B 645, 113 (2007). 

[7] E. Masso and J. Redondo, Phys. Rev. Lett. 97, 151802 

(2006) [arXiv:hep-ph/0606163]. 

[8] E. Zavattini et al. [PVLAS Collaboration], Phys. Rev. 
Lett. 96, 110406 (2006) [Erratum-ibid. 99, 129901 

(2007) ] [arXiv:hep-ex/0507107]. 

[9] See, J. Redondo, arXiv:0807.4329 [hep-ph], and refer- 
ences therein. 

[10] G. Cantatore [PVLAS Collaboration], Lect. Notes Phys. 
741, 157 (2008). 

[11] P. Sikivie, Phys. Rev. Lett. 51, 1415 (1983) [Erratum- 
ibid. 52, 695 (1984)]. 



[12] K. van Bibber, P. M. Mclntyre, D. E. Morris and 
G. G. Raffelt, Phys. Rev. D 39, 2089 (1989). 

[13] R. D. van der Hilst, et al, Science 315, 1813 (2007). 

[14] E. Clementi and D. L. Raimondi, J. Chem. Phys. 1963, 
38, 2686. 

[15] T. Nautiyal and S. Auluck, Phys. Rev. B 32, 6424 (1985). 
[16] G. G. Raffelt, Stars as Laboratories for Fundamental 

Physics (The University of Chicago Press, 1996). 
[17] G. G. Raffelt, Phys. Lett. B 166, 402 (1986). 
[18] H. Primakoff, Phys. Rev. 81, 899 (1951). 
[19] S. O. Costin and S. L. Butler, Phys. Earth Plan. Int. 157, 

55 (2006). 

[20] F. T. . Avignone, arXiv;0810.4917 [nucl-ex]. 

[21] K. Zioutas et al. [CAST Collaboration], Phys. Rev. Lett. 

94, 121301 (2005) [arXiv:hep-ex/0411033]. S. Andria- 

monje et al. [CAST Collaboration], JCAP 0704, 010 

(2007) [arXiv:hep-ex/0702006]. 
[22] See, for example: I. C. Atkinson et al., J. Magn. Reson. 

Imaging 26; 1222 (2007). 
[23] A. S. Chou et al. [GammeV (T-969) Collaboration], Phys. 

Rev. Lett. 100, 080402 (2008) [arXiv:0710.3783 [hep-ex]]. 
[24] A. J. Miller, et al, Appl. Phys. Lett. 83, 791 (2003). 
[25] A. E. Lita, A. J. Miller, and S. W. Nam, Optics Express 

16 3032-3040, (2008). 



4 



